Many-body dynamics of a Bose system with attractive interactions on a ring 
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We investigate the many-body dynamics of an effectively attractive one- dimensional Bose system 
confined in a toroidal trap. The mean-field theory predicts that a bright-soliton state will be formed 
when increasing the interparticle interaction over a critical point. The study of quantum many- 
body dynamics in this paper reveals that there is a modulation instability in a finite Bose system 
correspondingly. We show that Shannon entropy becomes irregular near and above the critical 
point due to quantum correlations. We also study the dynamical behavior of the instability by 
exploring the momentum distribution and the fringe visibility, which can be verified experimentally 
by releasing the trap. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in quasi-one- 
dimensional systems have attracted much more attention 
in recenWearsQ|l|l|l|l|![l3 

IIMIig . l2nll2lll2a |. because it, depending on the repulsive 
or attractive nature of two-body interaction, can form 
dark- or bright-soliton states, respectively 0, H, @- It 
is believed that atomic matter bright solitons are of pri- 
mary importance for developing concrete applications of 
BECs in the future Q, so of particular interests have 
been focused on the study of effectively attractive BECs 
II II Elll El E3 Using the Feshbach resonances [3, 
the coupling constant was tuned to a negative value and 
both bright-soliton trains ^3 an d a single bright-soliton 
|l5j were realized experimentally. 

The Gross-Pitaevskii mean-field theory (MFT) is a 
powerful tool in the study of BECs, which has suc- 
cessfully described an impressive set of experiments 
H3 IitI ITsf . However it breaks down if quantum many- 
body (QMB) correlations are important, such as in quan- 
tum phase transitions 0, H 0- For one-dimensional 
Bose system with attractive interaction, MFT predicts 
the transition from a uniform state to a soliton state [8| at 
a critical point when increasing the absolute value of the 
strength of the interaction. Recently, Kanamoto et.al. 
[l9| and Kavoulakis [2(j carefully checked the ground 
state and low-lying excitation spectra near and above 
the transition. They revealed that quantum correlations 
become crucial. The resulting properties of QMB correla- 
tions, such as the singularity in the transition is replaced 
by a crossover region, the condensate fraction begins to 
decrease, etc. [lji |2(|, suggest that the MFT should 
qualitatively be modified and QMB effects have to be 
considered. On the other hand, most of the previous 
studies regard QMB effects in the stationary situations, 
which attempt to beyond MFT by using the diagonaliza- 
tion schemes 0, |2(| • Therefore it is desirable to explore 
the quantum dynamics of this system. In the present pa- 



per, we numerically study the quantum dynamics of N 
interacting Bosons in a toroidal trap, which plays an im- 
portant role in the physics of trapped dilute gases |22{ . 
As a function of interaction strength, we calculate the 
Shannon entropy of the wave packets to characterize the 
quantum dynamics. A modulation instability at critical 
point is discovered, which is reflected by a rapid increase 
of entropy. As a direct result, the interference fringes, 
which occur after switching off the confinement and let- 
ting the particles spread in a ballistic way, vanish near 
and above the critical point. 

The toroidal one-dimensional regime can be realized 
when the transverse dimensions are on the order of the 
healing length, and the longitudinal dimension is much 
longer than the transverse ones H, 0, El El IH3 ■ Exper- 
imentally, this geometry may be achieved by an optical 
trap with Laguerre-Gaussian beams |2^|. The system is 
tightly confined in the radial direction so that the energy- 
level spacings in this direction greatly exceed the inter- 
action energy. Thus N Bosons are confined on a ring of 
radius R and cross section S = irr 2 with r <C R- The 
Hamiltonian for this system reads 0, Il9l | 



H 



U 



(1) 

where ^>(0) is the field operator, U = 8iraR/S, with a 
being the s-wave scattering length, and 9 denotes the az- 
imuthal angle. Here the length and the energy are mea- 
sured in units of R and ti 2 /2mR 2 . The normalization of 
*(6>) is J \ ^(6)\ 2 d6 — N. The Hamiltonian is character- 
ized by a single dimensionless parameter 0, El 
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which is essential the ratio between the interaction energy 
and the kinetic energy. The transition occurs when 7 < 

7c = -1/2 ii Mm 



2 



This paper is organized as follows. In sec. II, after 
numerically obtaining the eigenenergies and eigenstates 
of the system, we study the dynamical evolutions of 
Shannon entropy of the system. The evolutions show 
quasiperiodical and irregular behaviors below and above 
critical point correspondingly. The irregular one is inter- 
preted as much more number of virtually excited particles 
involved in the dynamics due to strong interaction. In 
sec. Ill, we discuss the momentum distribution of single- 
particle states. Above critical point the system exhibits 
a statistical relaxation to a steady distribution. The in- 
terference fringes of momentum, which can be obtained 
experimentally by releasing the trap, are studied too. We 
conclude in sec. IV 



II. SHANNON ENTROPY OF THE SYSTEM 




In cylindrical coordinate expanding the field operator 
^(0, t) in terms of plane waves as 

with I integer. Where c\ is annihilation operator of a 
boson with angular momentum I. Then Hamiltonian Q 
can be rewritten as 

H = Y, f2 ^ £l + 2N^ 6 i^fimC„ (4) 

l kl;mn 

Expressing the single-particle sates \l) according to 
their angular momentum I — 0,±1,±2,..., a 
many-body basis state \k) can be represented as 
, n_i, no, ni, " • • i n L,) restricted to the conser- 
vations of the total number of particles and the total 
angular momentum 

L L 

m = N,J2 Ith = (5) 

l = -L l=-L 

where ni denotes the number of bosons occupying the 
single-particle sates \l), and L is the cutoff of the angular 
momentum. Since N bosons define the smallest spacing 
on a ring, it is practical to make L fa N 22). The sys- 
tem is occupied by m — 2L + 1 single-particle states. 
With these assumptions, eigenvalues and eigenstates of 
Hamiltonian (0J can be solved directly by a smart method 
[24J (see Appendix). 

At first we introduce Shannon entropy to character- 
ize the dynamical behavior of the system, which has 
been used to manifest the many-body dynamics of Tonks- 
Giraudeau gas 22] . The Shannon entropy is defined 

S(t) = -J2\* k (t)\ 2 ln\* k (t)\ 2 (6) 

k 



FIG. 1: Entropy versus rescaled time t for L — 8, N = 8 
and different 7. Dotted curve describes 7 = —0.14. Dashed 
curve corresponds to 7 = —0.54. Thin solid curve represents 
7 = —1.62 and thick solid curve represents 7 = —6.48. 
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FIG. 2: Time dependence of principal components for same 
parameters as fig.l. The correspondences between line styles 
and 7 are same as fig.l too. 

Where ^(t) = (k\^(t)) is the projection of the wave 
function onto the many-body basis state \k). And the 
time-dependent wave function reads 

|*(t))=e-""|*(0)) (7) 

Where |^(0)) is state of the system at t = 0. 

Then, by varying a via the Feshbach technology [l4| . 
let us consider an interesting case for which initially the 
system is condensed on the ground state without interac- 
tion. After switch on the interaction, the wave function 
spreads over all the basis. 
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FIG. 3: Renormalized entropy for different 7 with N = 8 
and L = 6, 7, 8. The sharp increase indicates the instability 
region. Here we plot the absolute value of 7 for convenience. 



tropy is appropriate to characterize the dynamical prop- 
erties of the system. To compare the results for different 
7, we calculate the rescaled entropy S nor — S a ve/Njj, 
where S ave is the mean value of the entropy averaged 
over time after the saturation and Njj is the dimension 
of the Hilbert space, i.e. the total number of basis. Fig 
3 presents the results for L = 6, 7, 8 and N = 8. Here, 
different L can make us find the dependence of entropy 
on L. We can see that the data are sensitive to the modu- 
lation instability. For > 7 > y c , S nor is stagnated to 
zero. This is same to the MFT results, where bosons are 
almost condensed to the ground state. Near the critical 
point, the virtual high- momentum single-particle states 
are excited and the system is influenced significantly by 
the quantum fluctuations |19| . In this region, slight vari- 
ations of 7 alter the renormalized entropy obviously for 
given L. For 7 <C 7 C5 S nor reaches to a plateau with 
altitude ~ 0.85. We guess that for the limitation case 
L — > 00, the renormalized entropy will reach to unity 
corresponding to the situation that all of the basis states 
have the same probabilities to be occupied. 



The time dependence of entropy is reported in Fig.l. 
For > 7 > 7c, the entropy oscillates quasiperiodically 
in time. For 7 < j c , the small oscillation is disappeared 
and the entropy increases linearly to a saturation for a 
strong enough interaction. For weak interactions the sys- 
tem is condensed on the ground state. The small devi- 
ations from ground state are suppressed with time in- 
creasing 0, [2(| , which leads to the periodic oscillations 
of the entropy. It also ensures the validity of MFT with 
only a few modes [lj|,|2(j. When 7 is lowered (|7| is in- 
creased) , the number of virtually excited particles due to 
stronger interaction is increased 0, |2jj . This indicates 
that much more basis states are involved in the dynam- 
ics. As a result, the entropy is evolved with a linear 
increase followed by a saturation. It implies that the ini- 
tial ground state decreases exponentially while the basis 
excited in the dynamics increase exponentially in time 
(see Fig. 2). It can not be explained in the MFT frame- 
work. To address the number of basis obviously, one 
may calculate the principal components for wave pack- 
ets in the basis representation [25|. It is defined by the 
entropy N e ff(t) = exp[S(t)], which is regarded as the ef- 
fective number of unperturbed many-body states in the 
dynamics. A recent paper about the two-body random 
interaction model (23 shows that in a finite system, the 
number of basis states excited in the system is finite too. 
Although Hilbert space is same for the system indepen- 
dent of strength of interaction, the dynamics will depend 
upon it. Stronger the strength of interaction, more ba- 
sis involved in the dynamics. So the finity of basis leads 
to the saturation of the entropy. And the linear part of 
the time entropy may be interpreted as the onset of the 
chaotic superpositions of the basis states [2j| . 

Although the complicacy of the quantum correlations 
near and above the critical point [19j, the Shannon en- 



III. MOMENTUM DISTRIBUTION AND 
FRINGE VISIBILITY 

Quantum dynamics characterized by the Shannon en- 
tropy is investigated for different 7 in previous section. 
Similarly to MFT, a modulation instability is founded 
at a critical point. But near and above critical point, 
the irregular behaviors of quantum dynamics can not be 
obtained by MFT due to quantum correlations. 

QMB effects influence not only the entropy, but also 
the momentum distribution, which can be observed in 
experiment. What we do in the following is to address 
the influence. 

Using the formal Fourier transforming, we obtain the 
momentum function 

X(p) = J dri>(6)e-^ (8) 

Let us consider momentum along the y-axis p = p-e y . 
We get 

271 

x{p)=A^c k J d9e ikd e- ipRsine (9) 
k 

Where A is a normalization constant. Note that the 
integration is the Bessel function j2(| and substitute it in 
the equation. Then 

X (p) = A^CnMpR) (10) 

n 

Where J n (x) is the first kind Bessel function of order n. 
The occupation number distribution is 

n(p,t) = (y(t)\ X l(p)x(p)\y(t)) (11) 
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FIG. 4: Particle density distributions in the momentum rep- 
resentation for four parameters 7 at different times. Thicker 
curves describe latter times. Dotted, Dashed, Thin solid, 
thick solid curves correspond t=0, 0.1, 0.6, 1. 



Using Eq.J7J and Eq. ljTUj) , we obtain 

n(p,t) =Y J *l*vY; J ^ R ) J v(P R )( k \&vW) (I 2 ) 



k,k' 



LI' 



With (k\cjci'\k') — n^Skk'Sw where nf denotes the num- 
ber of particles with angular momentum I in the basis 
|fc), we get 




FIG. 5: Fringe visibility of occupation number for different 7 
. The curves correspond to same 7 as fig.l. Inset is the en- 
larging visibility for 7 = —0.14 , which shows clearly periodic 
behavior. 



condensates. This suggests that we can use the fringe vis- 
ibility to address dynamical evolutions. The visibility is 
defined as [22| 



v(t) 



(14) 



(13) 



Where 7 max = n(p = 0,t),I min = n(p 0: t) and p is the 
first zero of Jo- The results are reported in Fig. 5. The 
periodical and irregular behaviors of the dynamics are 
showed for some 7. For > 7 > 7 C , quantum dynam- 
ics study reveals that the visibility of a uniform ground 
state oscillates with time (see inset of Fig. 5), contrast- 
ing wi th the GP framework where the visibility is unity 
|a EH • For strong interaction, the visibility is described 
by a fast decay followed by irregular fluctuations. It can 
be used to distinguish the modulation instability. So en- 
tropy, momentum distribution and visibility are strictly 
consistent to characterize the dynamical behaviors of the 
one-dimensional Bose system. 



IV. CONCLUSION 



Where C is a normalization constant. 

The momentum distributions for several 7 are shown 
in Fig. 4. One can see that for weakly interacting system 
the occupation number distribution will not be changed 
obviously in time. But small perturbations lead to quasi- 
periodical oscillation of the entropy (Fig.l). When the 
two-body interaction strength exceeds the critical point, 
many high-momentum single-particle states are excited. 
With time increasing, the momentum distribution be- 
comes flat and the condensate can not be retrieval. A 
steady momentum distribution means the collapses of the 



We studied the dynamical properties of attractive one- 
dimensional Bose system with the QMB dynamics. Our 
main result is the onset of a modulation stability at j c . 
We found that the quantum dynamics of the system is 
distinguished by the periodical and irregular behaviors 
for 7 below and above critical point respectively. 

Within the MFT, a phase transition between uniform 
state and bright soliton is expected at a critical point || . 
We investigate the dynamical behaviors by varying the 
interacting strength. Far bellow the critical point, the 
system is confined to the ground state corresponding to 



a condensate. The entropy oscillates quasiperiodically in 
this situation, which is connected with the collapses and 
revivals of the condensate fraction (Fig.l and Fig. 2). And 
the system is well described by MFT. In the crossover re- 
gion the behaviors of entropy change gradually depending 
on 7. For strong interaction regime, evolution of the en- 
tropy is characterized by a fast linear increase followed 
by a saturation. This indicates that the ground state is 
relaxed to steady states distribution and the condensate 
is destroyed. The results of momentum distributions and 
the fringe visibility confirm the modulation instability of 
the system. They are of particular interest because they 
can be studied experimentally. 

The QMB effects, such as entropy, principal compo- 
nent, momentum distribution and visibility have impor- 
tant meaning in the study of the finite Bose system 
|22| . The investigation in the context demonstrates that 
QMB effects significantly influence the dynamical behav- 
iors and can help us to comprehend both the quantum 
evolution of the system and the mechanism of the tran- 
sition beyond MFT. 

Acknowledgments 

The work is supported in part by the NSF of China 
(Grant Nos. 60478029 and 10125419 ) and by the Na- 
tional Fundamental Research Program of China, Grant 
No.2001CB309310. W.L. thanks Dr. Y. Wu for a stimu- 
lating discussion. 

APPENDIX: ENERGY SPECTRUM AND 
EIGENSTATES OF A BOSE SYSTEM ON A RING 

In the following, we give an introduction about the al- 
gebra method to calculate the eigenvalues of Hamiltonian 
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(@J based on Ref. [24|. One can see Ref.0 for a detail. 

A many-body basis can be represented as 



\n-L,-" ,n_i,n ,ni, • • • ,n L ) = J— TT(cJ) n ' \ vac ) 

(A.l) 

where \vac) is the vacuum state satisfying ci\vac) = 0. 
Hence, we denote the eigenstate of Hamiltonian (@J) as 



I*} = F(cl L , • • • , cl v 4 4, • • • , c[)\vac) (A.2) 



where F is the linear combination of the terms Yii ■ 
The eigenvalue equation ff |\&) = E\^f) can be rewritten 
as EF\vac) = HF\vac) = [H, F]\vac) due to H\vac) = 0. 
Using [c u F] = OF/dcl [c/c m , F] = d 2 F/dc\dc* m , the 
polynomial function F satisfies the following differential 
equation 



/ kl;mn 

where X[ = cj. And the polynomial function F = 
F(x—i,,X—i,Xq,Xi,--- ,Xl) is the linear combinations of 
the terms Yli x i l ■ Now, the eigenvalues problem is re- 
duced to seek polynomial solutions to a second order 
liner differential equation. Utilizing this method, it is 
not a hard work to use a MATHEMATICA program to 
calculate all the eigenvalues and eigenstates of Hamilto- 
nian (0J for finite N and L. 
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